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Abstract. We model the present-day number and characteristics of Ultra-compact X-ray binaries (UCXBs) in the Galactic 
Bulge. The main objective is to compare the results with the forthcoming data from the Galactic Bulge Survey (GBS) in 
order to learn about formation scenarios, notably core collapse supernova rate, common envelope evolution and dynamical 
interactions, as well as the binary initial mass function. We use a binary population synthesis code and detailed stellar 
evolutionary tracks, combined with observations of the 11 known UCXBs. We predict a few hundred faint UCXBs with 
periods of 20 — 40 minutes that would be detectable by the GBS. 
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INTRODUCTION 

Ultra-compact X-ray binaries (UCXBs) are low mass X-ray binaries with an orbital period below ca. 90 min, indicating 
a compact and therefore hydrogen deficient donor [1]. Via Roche Lobe overflow (RLOF) the donor loses mass, which 
is partially accreted by a very compact companion (a neutron star or a black hole). Due to friction in the accretion disk, 
the gas heats up to temperatures corresponding to X-ray radiation. 

The Galactic Bulge is the best direction to look in to find UCXBs because of the high local star concentration and 
the fact that the line of sight from Earth crosses the Galactic Disk. Also, UCXBs are too rare for us to find nearby in 
the Milky Way, and too faint to resolve in other galaxies except for the closest. 

The Galactic Bulge Survey [2] (GBS) is an X-ray and optical survey focused on a 12 square degree region around 
the center of the Milky Way (fig. 1). The aim of this study is to predict and explain observations made by the GBS, 
and gain a better understanding of the physical processes in UCXBs. 
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FIGURE 1. The target area of the Galactic Bulge Survey consists of two strips of 6° x 1° around the Galactic Center. Image 
courtesy of P. Jonker. 



METHOD 



The binary population synthesis code SeBa [3] has been used to evolve a population of zero-age main sequence 
binaries. From the evolved population two classes of systems are selected, with donors being white dwarfs or heUum 
burning stars. Systems must have a neutron star or black hole as accretor. A third formation scenario involves main 
sequence donors and magnetic capture [4] but this class is excluded here due to the comparatively low number of 
ultra-compact systems that appear in this way. 

For the subsequent evolution more detailed stellar evolutionary tracks are used. For the white dwarf donors, a zero- 
temperature approximation [5] has been used after the star has cooled sufficiently. The tracks for helium buming donor 
are taken from [6]. 

The star formation history of the Galactic Bulge is subsequently used to arrive at the present-day number and 
characteristics of UCXBs in the Bulge. Using longterm observations by the Rossi XTE All Sky Monitor (ASM), these 
model parameters can be converted to observational parameters, which corresponds to a prediction of the present-day 
observable population. 

RESULTS 
Population synthesis 

Helium buming donor systems are typically formed less than 500 Myr after system formation, and start RLOF 
within another 200 Myr, while the majority of white dwarf donor systems start RLOF only after 1—6 Gyr. Part of the 
explanation is that helium buming stars evolve from secondaries with a zero-age main sequence mass of 3 — 4 Mq and 
therefore have a shorter main sequence timescale than the least massive white dwarf progenitors have. Also, if the orbit 
at the time of the formation of the helium burning star is too wide, it will turn into a white dwarf before starting RLOF, 
disquahfying itself from the sample and leaving only systems take reach RLOF relatively early. The primary explodes 
in a supemova early on, after typically a few x 10 Myr. Most of the white dwarf donor systems reach RLOF within 
2 Gyr, but a significant fraction of the systems take much longer This is the case for both helium and carbon-oxygen 
white dwarfs. Helium white dwarf donors typically arise from approximately 2 M© secondaries. Carbon-oxygen white 
dwarf donors mainly evolve from 3 — 10 Mq primaries or secondaries that do undergo a helium flash. UnUke heUum 
buming stars, white dwarfs can take arbitrarily long to reach RLOF. 

Many white dwarfs will merge upon RLOF (fig. 2). White dwarfs with a neutron star companion and a mass above 
0.83 Mq experience dynamically unstable RLOF, leading to runaway mass loss and a common envelope and eventuaUy 
a merger. For a black hole accretor, this Umit is approximately 1 .0 Mq . Of the remaining systems, the ones with a mass 
over 0.38 Mq cannot eject enough of the transfered matter from the system by isotropic re-emission. Even though the 
system is stable from the donor perspective, it is not from the accretor perspective and the system will still end up 
merging. In effect, almost all carbon-oxygen white dwarf systems merge, while almost all helium white dwarf systems 
do not merge. The exception are systems with black hole accretors. In that case low mass carbon-oxygen white dwarf 
systems can experience stable mass transfer. Systems that are unstable in either way are removed from the sample. 



Star formation history and birthrates 

Besides knowing which systems form and how long it takes for them to reach RLOF, we also need to know when 
these systems form: this is determined by the star formation history of the Bulge, shown in fig. 3. The star formation is 
mostly concentrated around 10 Gyr ago, but important is the tail of recent star formation in the C7 = 2.5 Gyr distribution. 

Combining the star formation history with the delay times of the start of RLOF yields birthrates for all systems, as 
shown in fig. 4. The main difference is that white dwarfs can start RLOF recently, while helium buming donor systems 
cannot, unless there is significant recent star formation (CJ = 2.5 Gyr). 

The evolution following RLOF or shortly before has been computed using tracks based on stellar evolution models 
such as [6]. For this specific stage of evolution, these tracks are more detailed than SeBa and therefore preferred. The 
evolution of UCXBs with white dwarf donors is governed by angular momentum loss by gravitational wave radiation. 
Mass loss from the donor is required to keep the radii of the donor and its Roche Lobe in accordance. After RLOF 
the donor radius will increase immediately for (fully degenerate) white dwarf donors and after roughly 100 Myr for 
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FIGURE 2. Total time-integrated amount of UCXBs with a white dwarf donor. Black: helium white dwarf donors. White: carbon- 
oxygen white dwarf donors. 




Time {Gyr) 



FIGURE 3. Star formation history of the Galactic Bulge as a Gaussian distribution for = — 10 Gyr and a = 
0.5 (solid), 1.0 (dashed) and 2.5 (dotted) Gyr. The total star formation is 1 • lO^" Mq. Note tiiat Time = corresponds to the 
present. [7, 8] 



helium burning donors (once the donor has become sufficiently degenerate following the extinction of nuclear fusion 
due to mass loss). The orbital period will increase with mass loss for degenerate donors. 

The amount of systems below a given period (fig. 5) initially follows the star formation rate as more systems are 
formed at short periods. As the star formation slows down, the amount keeps going up for a short while because 
more new systems are formed than older systems are removed from the given sample due to increasing periods. Given 
enough time, the amount will go down. Notably, the amounts of helium burning donor systems and white dwarf donor 
systems dechne at different rates corresponding to their recent birthrates (fig. 4). 



Present Bulge population 

Of special interest is the population today. For a star formation history distribution C7 = 1 Gyr, most systems are old 
and have expanded to a period of around 80 min (fig. 6). All short period systems are young and therefore must contain 
white dwarf donors, generally helium white dwarfs. White dwarfs with black hole companions reach periods longer 
than 90 min because the gravitational wave timescale is shorter for higher accretor mass, resulting in a higher mass 
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FIGURE 4. Birthrates for white dwarf donor systems (solid) and helium burning donor systems (dashed). Left: (7=1 Gyr. Right: 
a = 2.5 Gyr. 



FIGURE 5. Amount of UCXBs against time. Solid: white dwarf donors. Dashed: helium burning donors. = 1 Gyr. 

transfer rate and a lower donor mass than for a less massive accretor, after a given time. For o = 2.5 Gyr young helium 
burning donor systems exist as well, because of the recent star formation. Their donors have turned into carbon-oxygen 
white dwarfs during their evolution. 

We are interested in what we can observe and therefore must convert model parameters to observables. More 
specifically, we have to convert modeled mass transfer rates to observational luminosities. A reliable method is to 
use observations of known UCXBs by the ASM. The ASM has collected count rates in the 2 — 10 keV range since 
1996. These have been converted to bolometric luminosities using distance estimates from [9, 10] and a photon energy 
and bolometric conversion by [11]. Figure 7 shows how often a source emits above a given luminosity. Most UCXBs 
are highly variable. Most systems are only visible a small fraction of the time above a given luminosity, so at any 
given time we will only see a small fraction of the UCXB population using an instrument with a given sensitivity. The 
short period systems are generally visible more often than the long period systems, which is to be expected given their 
higher mass transfer rates. 

This data can be applied to the modeled systems. Depending on their period, a luminosity is assigned, drawn 
randomly from the time-luminosity curve (fig. 7) of the best-matching real UCXB. More than 99% of the modeled 
systems will be invisible at any given time. This causes a "population inversion": even though there really exist many 
more long period (around 80 min) UCXBs, from a observational point of view there are more short period systems. 
Applying the ASM data on the model results (fig. 8, a = 1 Gyr) gives the circles above the ASM sensitivity Umit. The 
data points represent the prediction of the present-day population of UCXBs. Some circles are below this limit because 
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FIGURE 6. Present-day orbital periods and mass transfer rates for white dwarf donors (black) and helium burning donors (white). 
Left: a = 1 Gyr. Right: a = 2.5 Gyr. 
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FIGURE 7. UCXB variability: fraction of time that a source emits at a given luminosity for 1 1 known UCXBs. [9, 10, 1 1] 



some real UCXBs are closer to us than the Galactic Center. In general though, below the ASM limit no prediction is 
possible using the ASM data. In a forthcoming paper these observational results will be explained in greater detail. A 
second method to convert mass transfer rate to luminosity is the Disk Instabihty Model (DIM) [11], which can be used 
to calculate luminosities fainter than ASM can observe. According to this theory, a mass transfer rate must exceed a 
critical value in order to be stable and the source to be persistent, i.e. visible at a high luminosity (almost) all the time. 
For lower mass transfer rates, the accretion disk builds until the cold outer disk becomes unstable and collapses onto 
the accretor. The source is only visible during these outburst stages, therefore the fraction of time the source is visible 
goes down for lower mass transfer rate, while the luminosity during these outbursts is held at a fixed value. 

At long periods, the actual number of systems will be much less than shown in fig. 8 due to tidal disruption. Tidal 
heating wiU change the mass-radius relation of the low mass donor leading to disruption of the system [12]. For 
a = 2.5 Gyr, there are more carbon-oxygen donors at 40 — 60 min which have evolved from heUum burning donors. 



DISCUSSION AND CONCLUSIONS 



• Tidal instabihty for low mass donors may be more effective than theory suggests, otherwise the predicted number 
of bright long period UCXBs is much larger than the observed. 

• Since the observed UCXBs (periods below 60 min) with known composition consist of helium and carbon-oxygen 
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FIGURE 8. Predicted observable population of UCXBs, based on ASM observations and DIM theory. C7 = 1 Gyr. 

in more or less equal numbers [13], recent star formation is required which contrains the star formation history 
model. A carbon-oxygen donor implies a helium burning donor origin because carbon-oxygen white dwarf donor 
systems merge upon RLOF in case of a neutron star accretor, and a black hole accretor system would likely be 
disrupted already unless it is very young. Short period (below 20 min) UCXBs are most likely helium white 
dwarfs because they have a higher formation rate at the present, regardless of which star formation history is 
used. 

• Prediction: the existence of a few hundred faint UCXBs with periods between 20 — 40 min that would be 
observable by the GBS. 
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